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' Abstract. The role of direct URCA reactions in damping of the gravitational radiation driven 

O^l i instability is discussed. The temperature at which bulk viscosity suppresses completely this in- 
stability is calculated. 

" ^ I The results are obtained analytically using recent calculations performed in the case of bulk 

\^ • viscosity due to the modified URCA processes (Lindblom 1995; Yoshida & Eriguchi 1995). 

I The bulk viscosity caused by direct URCA reactions is found to reduce significantly the region 

T-H ■ of temperatures and rotation frequencies where a neutron star is subject to GRR instability. 
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r-| . Introduction 

Oh| An absolute upper bound on angular velocity of rotating star is given by a shedding limit, ie. the 

Q ■ Keplerian frequency r^Kcpicr, defined as a velocity of a test particle on a circular orbit on equator. 

f~l , However it is possible that a star rotating with Q < fiKcpicr is unstable against some perturbations. In 

^ ■ this case this instability sets the upper limit on star's rotation. 

• • , Neutron stars rotating with very high angular velocity could be subject to secular instability caused 

' by the gravitational radiation reaction (GRR). This possibility has been pointed out by Chandrasekhar 

, (1970) and then shown to be generic feature of rotating star (Bardeen at al. 1977; Friedman & Schutz 1978; 

^ ; Friedman 1978). 

Further studies have shown that dissipation processes, such as bulk and shear viscosity, moderate this 
instability and in some ranges of temperature damp it completely. For temperatures lower than 10^ K shear 
viscosity (r/ T^^) stabilize stars rotating with < l^Kepier (Cutler & Lindblom 1987; Ipser & Lindblom 
1991). At high temperatures T > 10^° K the gravitational radiation driven instability is damped by bulk 
viscosity resulting from modified URCA processes (Cutler at al. 1990; Ipser & Lindblom 1991). 

It was recently pointed out by Lattimer et al. (1991) that the proton fraction inside neutron stars 
could be sufficiently large to make direct URCA processes possible. Haensel & Schaeffer (1992) calculated 
the bulk viscosity of hot neutron-star matter caused by direct URCA reactions. In the interesting range of 
temperatures this viscosity appears to be many orders of magnitude larger than bulk viscosity corresponding 
to modified URCA processes. Thus damping of GRR instability should be much more effective. 

In present paper we analyse the consequences of direct URCA reactions for the instability driven by 
gravitational radiation reaction. We use the results of previous calculations performed by Lindblom (1995) 
and Yoshida & Eriguchi (1995) for bulk viscosity due to the modified URCA processes. Our conclusions 
are obtained analytically using approximate formulae describing the damping of GRR instability for given 
value of viscosity. However, because of the strong dependence of shear and bulk viscosities on temperature, 
it seems that our estimates are very accurate. 
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Calculations and Results 

Bulk viscosity from direct and modiBed URCA processes. 

The bulk viscosity of the matter is given by the formula: 

C = ^-^^ (1) 

(Sawyer 1989) where n is the total baryon number density and eoefheient A describes the linear response 
of the rates of the beta and inverse-beta reactions to the departure from beta equilibrium (given by the 
value oi dfi — fin — ^J■p — ^J-c)- The coefficients B and C determine the change of Sfi due to the change in the 
chemical composition, at fixed density of the matter (B) and due to the compression or rarefaction of the 
matter, where the weak interactions are kept frozen (C). 

In our considerations we calculate bulk viscosity using "high-frequency" approximations both in the case 
of direct and modified URCA reactions. This approximations is valid if the second term in the denominator of 
Eq.(l) is negligibly small compared to w"^. The angular frequency of the pulsational modes which determine 
the maximum angular velocity of the star is of the order of 10* s~^. In our case we can estimate the value 
of (2AS/nu;)2 as: ~ 0.05 • T^^y/u;l for direct URCA reactions and ~ lO'^ • T^ly/tul for modified URCA 
processes (TMeV = kT/1 MeV and uj4 = w/lO'* s~^). We see that at, say, T < 10^° K we can safely use 
"high-frequency" approximation {{2\B / nco)"^ <t. 1). 

In this limit bulk viscosity caused by modified URCA processes is given by: 

U = 6.0 -101^(^)^96 gcm-is-\ (2) 

where pig = p/W^'' gem-^ and Tg = T/IO^ K. 

The above formula was used by Ipser & Lindblom (1991) and had been derived by Sawyer (1989) but 
published originally with a typographical error (Lindblom 1995). 

The bulk viscosity due to the direct URCA processes is given in the "high-frequency" limit by (Haensel 
& Schaeffer 1992): 

^^ = 8.9.10'''^^(x^Y'^T,' gcm-^s-S (3) 

where mn, m*, rUp, m* are neutron and proton true and effective masses respectively, Uq — saturation 
density of nuclear matter (no = 0.16 fm~^), x = rip/n — proton fraction, and (Cioo — in units 100 MeV). 
The parameter Cioo is density (and model) dependent and in the region of densities interesting for our 
problem is of the order of one (between 1 and 2, see Fig. 2 in the paper by Haensel & Schaeffer (1992)). It 
should be mentioned that to derive bulk viscosity of neutron-star matter Sawyer (1989) used free Fermi- 
gas (n-p-e) approximation and Haensel & Shaeffer (1992) take into account the effect of nueleon micleon 
interaction on the nuclear symmetry energy of the matter. These effects play a crucial role in determining 
the proton fraction x (sufficiently large value of x allows direct URCA reactions to take place) and also 
influence the function C{p). We disscus here two models corresponding to two different depencences of the 
nuclear symmetry energy on density (Model I and 11 of Haensel & Schaeffer (1992)). 
We will use a simplified version of formula (3) 

1/3 

Cd = 6.5-102Vof ^79* gcm-is-\ (4) 

where xqi = x/0.1 and we have introduced parameter q = ^^C'loc which is dimensionless and of the 
order of one. 

There are some physical limits on the region of densities and temperatures where direct URCA processes 
are operating and above formula is applicable. 

First, the value of the proton fraction x should be sufficiently large {x > Xcrit ~ 5)- This condition 
defines minimum density below which we have to do with modified URCA processes. This minimum density 
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Pmin was calculated by Hacnscl and Scliaoffcr (1992) and is equal to 0.53 • 10^''' g cm~"^ and 0.72 ■ 10"'^''' g cm^'' 
for their models II and I respectively. Thus direct URCA processes are operating in a central regions of a 
neutron stars. 

Second, there exist a critical temperature Tent below which neutrons become superfluid and protons 
superconductive. In this region reaction rates of the URCA processes (both direct and modified) are signif- 
icantly reduced. The estimation of Tcrit gives the value of the order 10® K (Page 1994). 

Damping of GRR instability by "direct URCA" bulk viscosity. 

To calculate the effect of bulk viscosity due to the direct URCA reactions on the damping of GRR 
instability we use previous results published by Lindblom (1995) and Yoshida & Eriguchi (1995) obtained 
taking into account modified URCA processes only. 

We restrict ourself to the Newtonian stellar model of the mass 1.5 Mq built of polytropic matter with 
index n = 1 (pressure P ^ p^). 

From physical point of view the mechanism of the damping of the oscilations which tends to grow 
as a reaction to gravitational radiation is in both cases (modified or direct URCA) the same. Namely 
the source of damping is departure from beta equilibrium due to compression or rarefaction of the matter 
during oscillations. The only difference are reactions between constituents of the matter which lead to above 
nonequilibrium processes. 

The conclusion is that for a given value of Cd the effect on damping of GRR instability is the same as 
for the same value of Cm- 

Above statement is true if we consider a small piece of matter in a star and would be true if we have 
constant viscosity C, throghout a star (or the dependence of Cm and Cd on p is the same). But this is not the 
case and damping of GRR instability by viscosity is a global feature of a star. The functional dependence 
C(p) enters ours considerations through the integral defining the damping time due to the bulk viscosity of 
matter (Ipser & Lindblom 1991): 

:^ = ^y'c5a,5a*(i-V (5) 

where a = Vqw" is the expansion (f" — fluid velocity). 

We take this fact into account defining the parameter C^: 

Q = ^/ i5a5a*d^x, (6) 

where Q = ({p) and p is the mean density of a star. It should be stressed that the value of depends on 
the pulsation mode. 

Using this definition we can rewrite Eq. (5) as 

i = Cc-C (7) 

Assuming, that for the same value of GRR instability is damped to the same extent we can very 
easily estimate the the temperature at which bulk viscosity caused by direct URCA reactions completely 
damps this instability (fimax = f^Kepier) having corresponding value Tm for modified URCA. The resulting 
formula is: 

logTd =1.5 log Tm - 6.26 + 0.25 log ( §^ 

5/3 (8) 

+ 0.25 log 

For modified URCA this temperature is equal Tm = 10^° K (Lindblom 1995; Yoshida & Eriguchi 1995). 
Thus in the case of direct URCA and polytropic star with index n = 1 and mass M = 1.5Mq considered by 
Lindblom (p^g = 0.36) we obtain: 
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logTd = 8.55 - 0.251og(gdm) - 0.25 log g^f (9) 

where we have introduced the parameter q^m = C'fd/^Cm- 

This parameter measures the relative difference between damping by modified and direct URCA reac- 
tions rcsiilting from a different shape of Cm and Cd throughout the star. 

To estimate the value of gdm we should evaluate two integrals: 

^•^ 2E J C 



Cd 

A more careful analysis of the radial dependence of integrands with the help of the Fig. 2 in the paper 
by Ipser & Lindblom (1991) leads to the following approximation of qdm- 

where 0d and are integrands in C^^ and C^^ normalised as in Ipser & Lindblom (1991) (the maximum 
value of each integrand is one) . 

The crucial point in estimating of C(^^ is the radius of the central sphere in the star where direct 
URCA processes are operating. Stricly speaking in calculation of C^^ we should restrict ourself just to this 
central region: 

^^^^^r^j"^"*'^'''' (12) 

or, in other words, put Cd = for r > fd. 

We can determine the value of rd for the considered model of a star (tt. — 1 polj^rope, /^central — 
1.2 ■ 10^^ gcm~'^). For models I and II of Haensel & Schaeffer (1992) we obtain: 

Model I pmin = 0.72 • lO^^ g cm'^ rai = 0.53 • R 

Model II pmin = 0.53 • IQi^ gcm-3 rdii = OM-R 
where R is the radius of the star, R = 12.53 km. 

Having the value of and comparing it with the Fig 2. of Ipser & Lindblom (1991) we could answer the 
question whether direct URCA processes are operating in the region where bulk viscosity damps effectively 
GRR instability . 

These two effects limit the region of a star where damping of GRR instability by direct URCA takes 

place. The oscillatory modes which are subject to this instability reveal significant changes in density in the 
outer region of a star. From the Fig. 2 in the paper by Ipser & Lindblom (1991) we could see that function 
5a Sa* starts growing at « 0.4 • R. 

Thus the main contribution to the integral C(;^ comes from the region between r^- and r^. 

This considerations enables us to estimate the value of qdm- <Zdm ~ 1/5-rl/lO for model II and 
9dm ~ 1/20 1/100 for model I. 

Taking into account these estimations we can rewrite Eq.(9) as 

log Td » 8.9 - 0.25 log qxli^ Model I 

log Td « 8.7 - 0.25 log qxli^ Model II 

We see that even though damping of GRR instability by direct URCA reactions is limited to rather 
narrow region of the star the process is effective for temperatures larger than 10^ K. 

Critical angular velocities of rotating NS. 

Using some approximations we can estimate not only the minimum temperature Td at which bulk 
viscosity caused by direct URCA processes completely damps GRR instability but also the function f2c(T) 
in the region where star is subject to this instability for sufficiently large fl {flc < O < fiKepier)- 



4 



For "modified URCA bulk viscosity" Oc(r) was determined by Liiidblom (1995) and Yoshida & Eriguchi 
(1995). In this case GRR instability is completely damped for temperatures lower than = 10^ K by shear 
viscosity and for temperatures greater than Tm = 10^'^ K by bulk viscosity. Thus for T <T^ and T > Tm we 
have ^^c{T) = riKepior- Between and Tm function ^c{T) is smaller than fiKepier having minimum equal to 
« 0.95f^Kepier at T^in « 2 • 10^ K (Fig. 1 of Lindblom (1995)). 

In the case of the direct URCA reactions ilc{T) differs from that of Lindblom (1995) in the region which 
corresponds to the damping due to bulk viscosity. This part of the function Qc{T) is increasing and reaches 
its constant and maximum value ^Kepler at temperature T = T^, significantly lower than T^. Furthermore 
the slope of flc{T)d is smaller than for modified URCA reactions flc{T)in- The reason is simply the different 
power in temperature dependence of Cd{T) and Cm(7')- From formula (8) we get: 

iT2) = :;i:jT-^] (Ti), (14) 



rflogTy^' 3 \dlogT 

and two derivatives are calculated for the same value of ilc (or for tq^(T2) = tq^(Ti)). 

Because of the rather strong dependence of ( and i] on temperature {(/t] ^T^) the region where shear 
and bulk viscosities play significant and comparable role is very narrow in temperature and so minimum of 
the function flc (T) is quite well defined by intersection of the parts of these dependences calculated for shear 
viscosity and bulk viscosity separately. 

More precisely we can determine the temperature Tmin at which ilc{T) has its minimum analysing the 
function: 

1 1 / 1 1 



where dimensionless functions /3(fi) and e(f2) were defined by Lindblom (1995) and Ipser & Lindblom (1991) 
and presented on graphs in their paper. Introducing, as in the case of bulk viscosity, the parameter C^: 

- = Cr,-V (16) 
Tr, 

we write the Eq.(15) in the form: 

m-{Cr,rj+~e{n)-CcO. (17) 



The estimates of parameters and C<; from Table 1 of Lindblom (1995) gives the ratio 7 = C^/C^ to be of 
the order 0.01 (7 depends on the mode and is equal: 0.67-10"^, 0.9 • 10~^, 1.2 • 10~^ for modes m = ^ = 4, 3, 2 
respectively). The difference between Cr, and reflects the fact that shear and bulk viscosities damps GRR 
instability in two different ways. For the same value of viscosity shear viscosity damps this instability about 
two orders of magnitude more effectively than bulk viscosity. 

We can make a remark that to estimate roughly the temperature where shear and bulk viscosity play 
comparable role one should find the crossover point of the functions r]{T) and 7 • C{T) rather than r]{T) and 

C(r). 

We see that to find the value of Tmin we have to analyse the function: 

f3{n)-(ri+~e{Q)-^0- (18) 
This approach leads to the following formulae: 

fog Tmin « 8.3 - i log qxli^ Model I 
fog Tmin « 8.2 - - fog qxy^ Model II 
To obtain the above result we have extracted the value of e(f2) from Fig. 13 of Ipser & Lindblom (1991). 
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The corresponding values of r^cmin could be estimated from Fig. 1 of Lindblom (1995): 



^cmin — 0.97 ^Kepl* 
^cmin = 0.975 r^Keplf 



ler 



ler 



Model I 



Model II 



(20) 



The above results seems to be accurate. This approximate method applied to bulk viscosity caused by 
modified URCA processes reproduces very well the value Tj^in obtained by Lindblom (1995). 

We conclude that if bulk viscosity due to the direct URCA processes is operating at temperatures greater 
than ~ 10® K the region in the T — Q plane where a star is unstable with respect to GRR instability is 
significantly smaller than in the case of modified URCA processes. 

Conclusions 

Our analysis shows a significant role of direct URCA processes in damping instability driven by grav- 
itational radiation of a rapidly rotating neutron stars. In the absence of other restrictions (superfluidity 
of neutrons, superconductivity of protons) "direct URCA" bulk viscosity is effective at temperatures larger 
than ~ 2 • 10^ K and completely damps GRR instability at ~ 5 • 10^ 10^ K. This process together with 
a lower bound defined by shear viscosity limits the allowed region where a rapidly rotating star could be 
unstable with respect to emission of gravitational waves to temperatures from the interval 10^ -r- 10^ K. To 
be subject to this instability a star should rotate with high angular velocity which appears to be very close 
to the Keplerian one — at most 2.5 3 % lower than the shedding limit. 

Taking into account the posibbility of superfluidity of neutrons and superconductivity of protons below 
critical temperature Tcrit would slightly change our conclusions. Theoretical calculations of Tcrit lead to a 
very uncertain and model dependent results (for a discussion see e.g., Page 1994). It seems that in the 
density range where direct URCA reactions are possible (p > Pmin) the critical temperature is determined 
by neutron pairing and lies between 10^ K and 6 • 10^ K for various theoretical models (Page 1994, and 
references therein). Superfluidity of neutrons and superconductivity of protons reduce URCA reaction rates, 
although it is not obvious whether bulk viscosity is switched off exponentially for T < TcHt (~ exp(— Tcrit/T")) 
or smoothly in a rather broad temperature region near Tcrit. Recent calculations performed by Haensel at 
al. (1995) support latter possibility and it seems that "direct URCA viscosity" could play a significant role 
even at temperatures few times lower than Tcrit (Haensel 1995). 

It should be mentioned that in the region where neutrons form a superfluid (at T < Tcrit) the shear 
viscosity of the matter results from electron-electron scattering and is of the same order as in the case of 
neutron-neutron collisions (Ipser & Lindblom 1991). Thus superfluidity of neutrons would change the shear 
viscosity much less than the bulk one and the net effect would be shifting of the temperature at which shear 
and bulk viscosity play comparable role to a higher value and lowering of the corresponding fJcmin- 

Although superfluidity of the matter supresses the role of viscosity in damping of GRR instability it 
also generates new dissipative mechanisms (for review see e.g. Pines & Alpar 1992). Recently Lindblom 
& Mendell (1995) using simple analytical model estimated the role of so called "mutual friction" , the most 
important effect in damping of GRR instability. They concluded that this dissipative effect completely 
suppresses the GRR instability below Tcrit- 

The next point which could influence our result is a size of the central region in a star where direct URCA 
processes are operative. On the contrary significant density changes in modes which becomes unstable with 
respect to gravitational radiation reaction are limited to the outer region of a star. Thus there exist a shell 
where direct URCA reactions could damp GRR instability. The thickness of this shell depends on the density 
profile througout a star (and so on equation of state) and also on the mass of the star. Our calculations 
has been performed for n = 1 polytropic star of a mass M = 1.5 Mq. For less massive stars density in the 
"damping region" could be so small that interesting shell would be very thin or even disappears. 

The analysis of the stability of a rotating star in the T — Q plane should take into account cooling 
timescale of this star. Strictly speaking a star becomes unstable with respect to GRR instability when the 
growing time of this instability iw short compared to the timescale at which a star would cool enough to 
leave the "instability region" (Yoshida & Eriguchi 1995). The cooling timescale is many orders of magnitude 
smaller for direct URCA processes than for the modified ones. Thus in the newly born neutron stars the 
interesting region will cool to 10^ K in minutes and to Tmin ~ 10^ K in days (Haensel & Schaeffer 1992). 
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Yoshida & Erigiidii (1995) published functions H„^ax{T) for several values of the growing timescale of GRR 
instability. Their results show that this effect would shift ricmin to a little higher value. 
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